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^^ In recent years, dynamic indentation tests have been shown to be useful both in 

identification of mechanical properties of biological tissues (such as articular carti- 
lage) and assessing their viability. We consider frictionless flat-ended and spherical 
I sinusoidally-driven indentation tests utilizing displacement-controlled loading pro- 

^ tocol. Articular cartilage tissue is modeled as a viscoelastic material with a time- 

OO independent Poisson's ratio. We study the dynamic indentation stiffness with the 

Q^ aim of formulating criteria for evaluation the quality of articular cartilage in order 

O to be able to discriminate its degenerative state. In particular, evaluating the dy- 

(T^ namic indentation stiffness at the turning point of the flat-ended indentation test, 

^^ we introduce the so-called incomplete storage modulus. Considering the time differ- 

^ ence between the time moments when the dynamic stiffness vanishes (contact force 

r* reaches its maximum) and the dynamic stiffness becomes infinite (indenter displace- 

• i-H ment reaches its maximum), we introduce the so-called incomplete loss angle. Anal- 

rN ogous quantities can be introduced in the spherical sinusoidally-driven indentation 

^ test, however, to account for the thickness effect, a special approach is required. We 

apply an asymptotic modeling approach for analyzing and interpreting the results 
of the dynamic spherical indentation test in terms of the geometrical parameter of 
the indenter and viscoelastic characteristics of the material. Some implications to 
non-destructive indentation diagnostics of cartilage degeneration are discussed. 

Key words: Viscoelastic contact problem, cartilage layer, dynamic indentation 
test, asymptotic model 
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Introduction 



Joint cartilage is known to have very limited repair capabilities and poorly regenerates. 
Intensive recent research and development have brought many innovations and also first 
clinical results in cartilage repair. However, recent clinical, radiological and histological 



evaluation techniques show somehow contradictory results (Kusano et al., 2011) and this 



is why measuring stiffness parameters of cartilage, especially in vivo measurements, are 
of novel interest nowadays. Cartilage stiflFness parameters can be measured in confined 
(Suh et al., 1995) or in unconfined (Armstrong et al., 1984) compression of cartilage 



specimen. However, both the confined and unconfined compression tests need sample 
preparation, usually cylindrically shaped specimens of cartilage, and therefore prohibit 
in vivo measurements. Furthermore, the mapping of the surface is limited by the sample 
size. These limitations are less restrictive than those usually encountered in indentation 
testing. 

The first mathematical model allowing to measure stiffness parameters of joint cartilage 
layer in indentation mode with fiat-ended as well as with spherical indenters was developed 



by Hayes et al. (1972|). In the case of a fiat-ended indenter of radius a pressed against a 



sample of thickness h^ the indentation stiffness defined as the ratio of the contact force P 
to the indenter displacement w is given by 



P 

w 



2aE 



1 



:l^c' 



(1) 



Note that compared to Hayes et al. (1972), we replace the shear modulus G with £'/(2(l + 
nu))^ where E is Young's modulus, u is Poisson's ratio. The Hayes model (1) is based 
on Hooke's law and takes into account the thickness effect through the stalling factor /^c- 
Because the widely used Hayes model assumes linear elasticity, it therefore does not take 
into consideration the fact that cartilage stiffness parameters are strain-rate dependent, 
and thus the Hayes model does not consider the dynamic nature of cartilage stiffness. 



Recall that for time-dependent materials (Tschoegl, 1997), the dynamic stiffness is charac 



terized by the complex dynamic modulus £"* = Ei + iE2 consisting of the storage modulus 
El and the loss modulus E2 with i being the imaginary unit (i^ = —1). On the complex 
plane, £"1, which is a real part of £"*, and £"2, which is an imaginary part of as imag- 
inary £*, represent the legs (catheti) of a right triangle with the hypotenuse of length 

l^*! = J El + E2. The loss angle, 5, results from the ratio of E2 and Ei through the 
relationship tan 5 = E2/E1. It should be emphasized that the storage and loss moduli Ei 
and E2 represent the response of a material to a sinusoidal loading scheme and actually 
depend on the corresponding angular frequency of sinusoidal oscillations u). In order to 
underline this fact we will write Ei{(jj) instead of Ei and so on. 

For in vivo (or ex vivo) measurements of cartilage stiffness and mapping a cartilage surface, 
a mechanical model for cartilage has to be prioritized considering dynamic properties of 
cartilage and measuring in indentation mode ("Ronken et al. , 2011). And, application of a 



single indentation test during arthroscopy allows one to evaluate the quality of cartilage 



and to detect osteoarthritic degenerative changes (Korhonen et al., 2003). The readily 



available systems on the market for arthroscopic measurements such as the 'Artscan 1000' 



(Lyyra et al., 1995; Toyras et al., 2001D allow only stiffness measurements of the cartilage 
surface without explicit considering Ei{(jj) and 6{(jj). 



To ascertain dynamic biomechanical properties of articular cartilage, Appleyard et al 



(2001) investigated a handheld indentation probe with a flat-ended cylindrical indenter 



operating in a vibration mode at a single- frequency of 20 Hz. Using the theory of |Hayes"et 
al. (1972), the absolute value of the eflFective complex dynamic modulus can be evaluated 



as follows: 



l^*HI = 



Po (1 - i^') 



(2) 



Wq 2a/^c 

Here, Pq is the contact force amplitude, Wq is the displacement amplitude. Note that 
the complex dynamic modulus £'*(cl;) is termed eflFective here, because in the case of a 
poroelastic material such as articular cartilage, the biomechanical response is dependent 
on the frequency uj and the boundary conditions for the sample as well. It was also 



observed (Appleyard et al., 2001) that when articular cartilage is indented at frequencies 



above 10 Hz there is marginal change in the eflFective parameters |£'*(cc;)| and 6{(jj) with 
the effective dynamic modulus being of similar magnitude to the 'instanteneous' elastic 
modulus generated during a rapid load step indentation test. Nevertheless, to the best of 
our knowledge, no study has analyzed thus far the relationship between the parameters 
of time-dependent materials measured in a vibration indentation test and in a single 
indentation test. In order to facilitate such a comparison, we consider sinusoidally-driven 
displacement-controlled indentation tests. Note that as a flrst approximation, the half- 
sinusoidal indentation history can be used for modeling impact tests. 

Measuring stiffness parameters of cartilage in indentation mode with spherical tipped 
indenters has advantages as well as drawbacks. At one hand, with spherical indenters the 
error obtained when hitting the surface not exactly perpendicular is much smaller than 
with flat-ended indenters, where the surface is touched with one edge of the indenter flrst. 
For example, when the surface with a spherical indenter will be hit with 80° instead of 
90°, the result for the stiffness will be underestimated by less than 2%. On the other hand. 



spherical indenters underestimate the inhomogeneity and changes in stiffness (Schinagl et 
al., 1997) as function of the indentation depth. This is why both theories, for flat-ended 



and for spherical tipped indenters are provided. 

As measuring mechanical properties gained a new importance in recent years, because dy- 
namic indentation tests have been shown to be helpful both in identiflcation of mechanical 



properties of articular cartilage and assessing its viability (Bae et al., 2003; Broom and 



Flachsmann, 2003). Indentation stiffness is now accepted as a fundamental indicator of 



the functional mechanical properties of articular cartilage ( |de Freitas et al. , 2006). The 
dynamic stiffness is deflned as the ratio of input force, P(t), to output displacement, w{t). 
Thus, for a time-dependent material like articular cartilage, the indentation stiffness de- 
pends on the indentation protocol, and generally it is a function of time. It is also well 
known that the indentation stiffness depends on the indenter size as well as on the sample 
dimensions (see, e.g., Eq. (1)). This follows from a comparison of the stiffness dimension 
MT~^ with the dimension ML~^T~^ of Young's modulus. From a geometrical point of 
view, articular cartilage is usually considered as a layer of constant thickness, h. In view 
of the relative mechanical properties of cartilage and subchondral bone, it is assumed that 
the layer is flrmly attached to a non-deformable base. Thus, the indentation scaling factor 
will depend on the aspect ratio a = a/h^ where a is the radius if the contact area. 



The above simple analysis is applicable for the linear relationship between the contact force 
P{t) and the indenter displacement w{t)^ where the contact radius remains unchanged in 
time. In spherical indentation, the force-displacement relationship requires a more acute 
analysis. It will be shown that the results of the dynamic spherical indentation of a time- 
dependent material depend on the level of indentation. 



To a first approximation (Hayes and Mockros, 1971; Parsons and Black 



1977; Lau et 



al., 2008), cartilage tissue can be evaluated mechanically as a viscoelastic material with 



a time-independent Poisson's ratio, z/, such that the overall constitutive behavior of the 
material is expressed in terms of its complex modulus £'*(cl;). Indentations tests for vis- 



coelastic materials were studied in a number of publications (Oyen, 2005; Cheng and Yang 



2009; Argatov and Mishuris, 2011j). We consider fiat-ended and spherical indentation tests 
utilizing displacement-controlled loading protocol with the indenter displacement modu- 
lated according to a sinusoidal law at an angular frequency u) = 27rf (rad/s), where / 
(Hz) is the loading frequency. We apply an asymptotic modeling approach for analyzing 
and interpreting the results of the dynamic spherical indentation test in terms of the ge- 
ometrical parameter of the indenter (indenter radius, R) and viscoelastic characteristics 
of the material. In particular, we examine the relationships between the storage modulus 
Ei{(ju) and loss angle S{(ju) and the so-called modified storage modulus E^^^i^^ ^o) and the 
modified loss angle S^^2{^) i^ the displacement-controlled sinusoidally-driven indentation 
test. 



The rest of the paper is organized as follows. In Section 1, we consider cylindrical fric- 
tionless indentation of a viscoelastic layer. In particular, the linear force-displacement 
relationship is outlined in Section 1.1, while the indentation scaling factor for the cylin- 
drical indenter is considered in Section 1.2. Based on the analogy with the case of harmonic 
vibrations (considered in Sections 1.3 and 1.4), in 1.5, we introduce the incomplete storage 
modulus and loss angle as material characteristics that can be assessed directly from a 
single sinusoidally-driven indentation test. 

In Section 2, we study spherical indentation of a viscoelastic layer. Based on the general so- 



lution obtained by Ting ( 1968), in Sections 2.1 and 2.2, we write out the force-displacement 



relationship for the loading and unloading stages, respectively. The indentation scaling 
factor for the spherical indenter is introduced in Section 2.3. In Section 2.4, we introduce 
the so-called modified incomplete storage modulus and loss angle, and investigate their 
behavior in Section 2.5 for the standard viscoelastic solid model. 

In Section 3, we actually consider the thickness effect in spherical indentation of a vis- 
coelastic layer. By analogy with the elastic case, we introduce the quantity Ey^i^^^o) 
(which is called the modified storage modulus, in view of its relation to the storage mod- 
ulus Ei{(jj)) while the modified loss angle ^3/2(^) is introduced according to a standard 
interpretation of the time lag between the peak force and peak displacement. Some prop- 
erties of these parameters, which turned out to be dependent on the level of indentation, 
are illustrated for the standard viscoelastic solid model. Low- and high-frequency asymp- 
totic analysis of the quantities E^{(jJ^zuo) and S^{uj) is presented in Sections 3.2 and 3.3, 
respectively. 

Finally, in Sections 4 and 5, we outline a discussion of the results obtained and formulate 
our conclusions. 



1 Cylindrical frictionless indentation of a viscoelastic layer 



lA Liner force- displacement relationship 



We consider a viscoelastic layer bonded to a rigid substrate indented by a flat-ended cylin- 
drical indenter. For the sake of simplicity we neglect friction and assume that Poisson's 
ratio, z/, of the layer material is time independent. Then, applying the elastic-viscoelastic 



correspondence principle ( [Christensen , 1971), one can arrive at the following equation 
between the applied force P{t) and the displacement of the indenter w{t) (Zhang and 



Zhang, 2004; Cao et al, 2010): 



P{t) = Y^^ >'c{<^) J E{t - t) — {t) dr. 



(3) 



Here, a is the radius of the contact area, h is the layer thickness, E[t) is the relaxation 
modulus, t is the time variable, t = 0— is the time moment just preceding the initial 
moment of contact, /^c(<^) is a dimensionless factor, which is determined from the solution 
of elastic contact problem for a cylindrical indenter, a is the relative radius of the contact 
area that is, i. e.. 



Note that the dependence of /^c(<^) on Poisson's ratio is not indicated explicitly Fig. la 
illustrates the behavior of /^c(<^) for diflFerent values of u based on the numerical solution 
obtained by Hayes et al. ( 19721 ). 
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Figure 1. (a) Indentation scaling factor for the cylindrical indenter as a function of relative 
contact radius; (b) Relative error of of the asymptotic approximation (7). 

Denoting E{t) = £'oo^(t), where £"00 is the relaxed elastic modulus (the limit of modulus 
E{t) at t ^ oc), ^(t) is the relaxation function, we rewrite Eq. (3) in the following form: 



dw 






(5) 



Inverting the relationship (5), we obtain 

w{t) 



2aEoo /^c(<^) 



1 f ^, ^dP , . ^ 



(6) 



where $(t) is the creep function. 



1.2 Indentation scaling factor for the cylindrical indenter 



According to Vorovich et al. (1974); Argatov (2002), the foUowing asymptotic model takes 



place for the indentation scaling factor ^d^)'. 
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Here, Oq and Oi are asymptotic constants depending on Poisson's ratio v given by 

(-ly 



(7) 



22^(m! 



/[l-/:(A) 



A'^dA. 



In the case of a layer bonded to a rigid base, we have 

2xsh2A-4A 
^ ^ " 2xch2A + l + x2 + 4A2' 

where x = 3 — 4z/ is Kolosov's constant. 

To determine the range of validity of the asymptotic model (7), we compare its predictions 
with the numerical solution given by Hayes et al. ( 19721 ). As it could be expected (see 
Fig. lb), the accuracy of the approximation i<i^^[a) given by (7) decreases as Poisson's 
ratio approaches 0.5 . Fig. lb shows that asymptotic approximation (7) is quite accurate 
in the range a G (0,0.6), that is for the indenter diameter less than the layer thickness. 
Note here that the substrate effect on the incremental indentation stiffness was considered 



in the elastic case in (Argatov, 2010) 



1.3 Harmonic vibration 



Observe that Eq. (3) assumes that the layer material was at rest for t < 0. In order 
to study harmonic vibrations of a viscoelastic layer, we should replace Eq. (3) with the 
following one: 

^/ X 2o ,,/"_, .dw,^^ 

Pit) = ^^^^'^cl^) J E{t-T) — {T)dT. (8) 



Substituting a harmonic displacement w{t) = Im{t(;oexp(icL;t)} with amphtude wq and 
frequency uj into Eq. (8), one can arrive at the following equation: 

Pit) = Y^ '^'^^''^ lTn{E*{u) expM)}. (9) 

Here, Im denotes the imaginary part of a complex number, E*{(jj) is the complex relaxation 
modulus given by 

oo 

£;*(cj) = iuj f E{s)exp{-iujs)ds. (10) 



By convention ( Pipkin[[l986 ; Tschoegl, 1997), we define the storage modulus, Ei{laj)^ and 



the loss modulus, £2(00)^ as the real and imaginary parts of E^'^uj)^ respectively, i.e., 

E%uj) = Er{uj)+iE2{uj). (11) 

From Eqs. (10) and (11), it follows that 

00 
Ei{(jj) =(jjEoo hif{s) sinus ds, (12) 



00 
E2{u)) =ujE^ / "^^{3) C08 ujsds. (13) 



Furthermore, according to Eq. (9), we can write 

P{t) = Posin(cjt + 5), (14) 

where Pq is the force amplitude, 6 is the phase angle between the harmonic displacement 
and the force, given by the formulas 

^"^ K<,{a)\E*{u)\wo, (15) 



1-u 



Observe that the phase angle 6 depends on the frequency uj (this is not indicated in 
notation for simplicity). 

Finally, note that the vibration indentation tests should be accomplished with a quasistatic 
preload to ensure a complete contact between the indenter's base and the layer surface, 
since tensile stresses are not allowed in frictionless indentation. 



1.4 Determination of the complex relaxation modulus via vibration indentation tests 



We assume that the displacement and force amplitudes wq and Pq as well as the phase 
angle 6 are experimentally measurable quantities. Then, Eqs. (15) and (16) yield the 



following equations (Cao et al., 2010) 
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2a/^c(<^) '^0 



cos 5, 



sin 5. 



(17) 
(18) 



Thus, for a given constant frequency cj, the vibration indentation test yields the storage 
and loss moduli Ei{uj) and E2{(jo)^ if the amplitude ratio Pq/wq and the phase angle 5 are 
known from the experiment. 

Further, let t^ denote the moment of time when the indentation speed w{t) vanishes, 
that is, when w{tm) = and t^ = 7r/{2uj) + nk/u)^ /c = 0, ±1, ±2, .... Considering the 
indentation process over a half of period t G (0,7r/cL;), we will have t^ = 7t/{2(jj) and, 
correspondingly, w{tm) = '^o and P{tm) = PqCos^. Hence, taking Eq. (17) into account, 

we obtain the formula 

1 - z/2 P(t^) 



2aK^[a) w(tm) 
where t^ is a time moment such that w{tm) = 0. 



(19) 
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Figure 2. (a) Displacement-controlled oscillation test; (b) Displacement-controlled indentation 
test. 

Thus, according to Eq. (19), the ratio P{tm)/w{tm) at the time moment of the displace- 
ment extremum determines the storage modulus. 

Let now t^ be the moment of time when the indentation displacement w{t) vanishes, 
i. e., w{t'^) = and t^ = nk/uj^ fc = 0, ±1, ±2, .... (Note that for harmonic vibrations 
tm = tm + 7r/{2uj).) Taking into account Eq. (18), we obtain 



Eo(uj) 



1 - V^ UjP{f^ 



2a/^c(^) ^(C) ' 
where w{t'^) is the indentation speed when the indentation displacement vanishes. 



(20) 



Remark 1 For the sake of completeness, we provide below the dual- conjugate formulas 
for Eqs. (19) and (20). Let tM denote the moment of time when the derivative of the 



contact force, P{t), vanishes, that is, when P{tM) = and tM = (7r/2 — 6)/(jJ + Trk/uj, 
/c = 0, ±1, ±2, .... Let also t'^ he the time moment when the contact force P{t) vanishes, 
i. e., P{t'j^) = and t'j^ = (tt — 8)/uj + nk/uj, A; = 0, ±1, ±2, .... Then, according to 
Eqs. (17) and (18)^ the following relationships hold true: 



Ei{uj) 2aK(^{a) w{t 



M) 






(21) 



(22) 



\E*iuj)\^ 1-Z.2 p(t'^) 

Recall that the magnitude of the complex modulus, |£'*(cl;)|^ is determined by the last 
formula (16). 



1.5 Indentation test with a sinusoidal displacement. Incomplete storage modulus and loss 
angle 



Let us first consider a single indentation test with a prescribed sinusoidal displacement 
according to the law 

w{t) = Wosinujt^ t G (0,7r/6L;). (23) 

Here, Wq is the maximum depth of indentation, cj is a given quantity having the dimension 
of reciprocal time. The quantity 

has a physical meaning of the time moment when the indentation displacement reaches 
its maximum. We emphasize that due to viscoelastic properties of the layer material, the 
duration of contact will be less than tt/uj. 

According to Eq. (3), we get 

2a y 

P{tm) = _ 2 f^c{c^)^Wo / E{tm " t) COS cjt rfr, (25) 



where tm = 7r/{2(ju). 

By analogy with Eq. (19), we define 

^i(^) = 7^ r^^ — f 5 (26) 

where w{tm) = '^o (see Eq. (23)). 

It is clear that the quantity Ei{(jj)^ introduced for single indentation test, differs from the 
storage modulus Ei{uj)^ introduced for vibration indentation test. 

In view of (25), Eq. (26) yields 

7r/(2cj) 

Ei{uj) = uj I E{s)8mujsds. (27) 



Recalling the notation E{t) = £"00^ (t), we rewrite Eq. (27) in the form 

7t/{2uj) 

Ei{cju) = uuEoo / ^(s)sincL;5(i5. (28) 



Comparing Eqs. (12) and (28), we see that their right-hand sides differ only by the integral 
upper limits. 

Now, let tM be the time moment when the contact force (3) corresponding to the inden- 
tation law (23) reaches its maximum, i. e., P(tM) = 0. Then, by analogy with the case of 
linear harmonic vibrations (see Eq. (14)), we put 

5(cj) = |-c^iM. (29) 

The quantity 6{(jj) is called the incomplete loss angle determined from the sinusoidally- 



driven displacement-controlled cylindrical indentation test (Argatov, 2012). 
In view of (24), formula (29) can be rewritten as follows: 

5(^) = 1 (^"^ - ^^) . (30) 

Z 6777. 



The interrelations between the quantities Ei{uj)^ 5{uj) and Ei^u))^ S{uj) were investigated 
in (Argatov, 2012^ . It was shown that they asymptotically coincide, respectively, in both 



the low and high frequency limits, while within the intermediate range of cj, the differences 
depend on the viscoelastic model in question, that is on the properties of the relaxation 
modulus E{t). 



2 Spherical frictionless indentation of a viscoelastic layer 

2.1 Force- displacement relationship in the loading stage 

Applying the general solution obtained by jTing (1968) for a class of viscoelastic contact 



problems in terns of the corresponding elastic solutions, we will have 

^^'^ = 3(1 -Xr I ^{'-irr^Hr))}n - r) dr, (31) 

^ ' 

Ht) = ^^'^QiMt)). (32) 

Here, a{t) is the variable relative radius of the contact area, i.e. (cf. Eq. (4)) 

a{t) = ^, (33) 



while J^{a) and Q{oi) are depending on Poisson's ratio v dimensionless factors such that 

J'(O) = 6^(0) = 1. 

We wiU use formulas (31) and (32) under the assumption that the relative contact radius 
a{t) monotonically increases in the time interval (0,t^), where t^ is a certain moment of 
time. 



2.2 Force- displacement relationship in the unloading stage 



Let us assume that the relative contact radius a{t) decreases to zero in the interval 
t G {tm-itc)-) where tc is the time of contact of the indenter with the layer surface. Once 



again, making use of the general solution derived by Ting (1968) for the case when the 



variation of contact radius posses a single maximum, we obtain 

Pit) = ^^^3^ J J^{c^ir?nMr))}^{t - r) dr, (34) 



w{t) = ^L{t)'g{a{t)) 

-/*(t-r)|: / ^{a{vrGHv))}^{r-v)dvdr\. (35) 

tm tiir) ) 

Here, ti(r) is the time moment prior to t^ such that the contact radius a{T) is equal to 
the prior contact radius a(ti(r)). The function ti{t) remains to be calculated. 

If the indenter displacement w{t) is prescribed, taking into account the relation a{ti) = 
a{t) for ti < t^ < t, we arrive at the equation 

w{t) = ^^^^g{a{t)), (36) 

from which we get 

ti(t)=^"'(^^^^^(o^(t))), tG [t^,te]. (37) 



Finally, note that in the inner integral in (35), for rj such that ti(r) < rj <tm^ the function 
a{r]) should be calculated according to Eqs. (31) and (32). 



2.3 Indentation scaling factor for the spherical indenter 

In the elastic case, according to the notation used in Eqs. (31) and (32), we have 

a'T{a), (38) 



3(l-i/2)i? 

w = —a^g{a). (39) 

R 

Here, a is the relative radius of the contact area as defined by formula (4). 
Representing Eq. (39) in the form 



^/wR 



aJg{a), (40) 



h 
and taking into account that ^(0) = 1, we see that Eq. (39) can be inverted as 



VwR , . .^_ 

a=^—g[w), (41) 

where 

\/wR 

- = nr- <*2' 

Substituting the expression (41) into Eq. (38), we obtain the following relationship: 

Here we introduced the notation 

j{w)=g{wfT{wg{w)). (44) 

It is clear that /(O) = ^((0) = 1. 
Finally, in view of Eq. (40) , we can represent Eq. (43) as 

P=^^^fl'Kla\ (45) 

where we introduced the indentation scaling factor 

K,{a) = f(a^Q{^)). (46) 

We emphasize that k^{ol) is normalized in such a way that /^s(0) = 1. Numerical values for 
K^ioi) for a range of parameters a and v are given in Table 1 based on the results obtained 



by Hayes et al. (1972). We note that k^{ol) = (3/2) y^/^, where x ^^d k are parameters 



employed in their analysis. 



According to Argatov (2001), the following asymptotic expansion holds true: 



Table 1 

Values of Ks{a) for the spherical indenter. 
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(47) 



Now, employing the forth-order asymptotic model constructed by Vorovich et al. (1974); 



Argatov (2002), we obtain 



K^[a) = 1 + a h a —^ 



a 



TT 



TT^ 



— a 



22ao 32aoai 



Stt^ 



157r2 



16al 32ai 



(48) 
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Figure 3. (a) Indentation scaling factor for the spherical indenter as a function of relative contact 
radius; (b) Relative error of of the asymptotic approximation (48). 

Fig. 3a shows details of the behavior of ^^{a) for different values of Poisson'a ratio. The 
errors of the asymptotic approximation K^^{a) given by (48) are plotted as functions of 
a/h in Fig. 3b based on the data given in Table 1. 



2.4 Modified incomplete storage modulus and loss angle 



In the viscoelastic case, according to Eqs. (31), (32), (42), and (43), we will have 



p{t) 



AE^h^ r d 



3(1 - iy^)R 



lj^{zu{Tf/'f{zu{T))}^{t-T)dT. 



(49) 



Here we used the notation (cf. (42)) 



zu{t) 



w{t)R 
It ■ 



(50) 



Equation (49) can be simplified for the case of a viscoelastic half-space when /(tu(r)) = 1 
as follows: 

4.^/R } d 



^("-3(f^/|{"M"^}^"-^'*- 



(51) 



We emphasize that Eqs. (49) and (51) are valid under the assumption that the indenter's 
displacement w{t) increases in the time interval (0,t^). 

Further, we consider again the same single indentation test with a prescribed sinusoidal 
displacement according to the indentation protocol (23). We may use Eqs. (49) and (51) 
in the time interval (0, tm) with tm = 7t/{2(jj)^ that is up to the moment, when the indenter 
reaches its maximum indentation depth Wq. 



By analogy with Eq. (26), we consider the quantity 



3(1 - Z.2) p(t^) 



4v^ w'J' , 



J E{t - t) 



d /w(r)\3/2 



dr. 



(52) 



dr V W() 
Note that the quantity on the right-hand side of Eq. (52) was previously considered in 



a number of studies on indentation of viscoelastic materials (Hu et al., 2001; Kren and 



Naumov, 2010). 



Substituting the expression (23) into the right-hand side of Eq. (52), we arrive at the 
following integral with /3 = 3/2: 



7r/(2cj) 



^'^('^>= / ^(i-^)5;(^'"^^''''^- 



(53) 



Observe that here the parameter /5 was introduced to simplify notation. However, later 
we show (see Remark 2) that the notation Ef^{uj) is meaningful for different values of /5. 



By changing the integration variable, the integral (53) may be cast in the form 

7r/(2cj) 



E^{uj)=uj I E{s)f3cos^ ^ ujssinujsds. 





(54) 



It is clear that for ;5 = 1, the right-hand sides of (27) and (54) coincide. The quantity 
Ei^{uj) will be called the modified incomplete storage modulus. 



Remark 2 Recall (Galin, 1946; Borodich and Keer, 2004) that the force- displacement 
relationship in the elastic case for a rigid blunt indenter with the shape function z = Ar^ 
is given by the equation P = [E/{1 — v'^)\A^~^Kf^w^ with /3 = (A + 1)/A and (with T[x) 
being the Gamma function) 



K,= 



22(^-i)(/5-l) 



-r 



,exp,l±Az^..(.-l)^-^^<^-': 



2(1-/5) / 1 x/5-1 



For a spherical indenter^ we have A = 2 and /3 = 3/2. We refer to (Argatov, 2011) 
for complete details of this consideration in the elastic case. In the viscoelastic case, the 
force- displacement relationship in the loading phase is given by 



Pit) 



EA^-^Kh r d 



1 



j j-{w{Tf]E{t-T)dT. 



Comparing this equation with Eq. (51), we see that the blunt indentation test yields the 
modified incomplete storage modulus Eis{u) introduced by formula (53). This explains the 
introduced notation. 



Further, assuming the variation of the indenter displacement in the form (23), we get the 



following variation of the contact force: 



4:VR 3/2 r ^., . d 



^ ^ 

Now, replacing 3/2 with /5 in the exponent under the integral sign in (55), we obtain 

P{t) = -^^^^wl'^ (3uj f E{s){smuj{t - s)f-^ cosuj{t - s) ds. (56) 

3(1 — u"^) J 



Now, let tM be the time moment when the contact force (56) reaches its maximum, i. e., 
P{tM) = 0. Then, by analogy with the case of linear harmonic vibrations, we put 

^/3(^) = 2 ~^*M. (57) 

The quantity Sis^u)) will be called the modified incomplete loss angle determined from the 
sinusoidally-driven displacement-controlled spherical indentation test. In view of (24), 
formula (57) can be rewritten in the form (30). 



2.5 Modified incomplete storage modulus and loss angle. Standard viscoelastic solid model 



In order to fix our ideas, we assume that the layer's material follows a standard linear vis- 
coelastic solid model, which is described by the following normalized creep and relaxation 
functions: 

$(t) = 1 - (1 - p) exp(-t/r,), ^(t) = !-(!- 1/p) exp(-t/(pr,)). (58) 

Here, Tg is the characteristic retardation or creep time of strain under applied step of 
stress, p is the ratio of £"00 to the unrelaxed elastic modulus Eq (modulus E{t) at t = 0), 
i.e., p = Eoc/Eo < 1. 



The following relations are well known (Tschoegl, 1997): 



Ei(c.) = E^ + {E, - E^) ^^if^, , (59) 

E2{uj) = {Eo-E^) 



u^pTsY + r 

5{u) = arctan^^-Jp^^. (60) 

Fig. 4a shows the behavior of the modified incomplete storage modulus Ej^iu) in com- 
parison with that of the storage modulus Ei{uj) given by (59) for /3 = 1, 1.5, and 2. At 
high frequencies, the both dimensionless quantities Ei{(jj)/Eoq and Ei^{uj)/Eoq approach 
the limit value E^/ E^ = 1/p = 2. Fig. 4b shows the relative error of the approximation 
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Figure 4. Modified incomplete storage modulus E^{(jj) 
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of Ei{uj) by £"3/2 (^) for different values of the dimensionless parameter p. In each case, it 
is assumed that the mean relaxation time pr^ is the same. 



Fig. 5a presents the behavior of the modified incomplete loss angle 5m (^) determined 
from the displacement-controlled indentation test in comparison with that of the loss 
angle S{(jj) given by (60) for /5 = 1, 1.5, and 2. The error of the approximation of the loss 
angle S{uj) by 5m (^) is shown in Fig. 5b. Unfortunately, the deference between 5m (^) and 
6{(jj) does not vanish as cj ^ oc. It will be shown that the limit value of the relative error 
is 33.33%. 



3 Accounting for the thickness effect in spherical indentation of a viscoelastic 
layer 



3. 1 Dynamic parameters for assessing the mechanical properties and viability of articular 
cartilage by a spherical indentation test 



Finally, let us consider the quantity 

3(1 - u') P{tJ _ ^0 

where /^s(<^m) is the indentation scaling factor corresponding to the maximum indentation 
depth am = Ci{tm)^ thus am = ctm/h. 



According to Eqs. (46) and (49), the right-hand side of Eq. (61) is determined as follows: 

^^{uj,wo) = jv^ J ^(^-rV{(sinc^r)^/V(^oVsinc^r)}rfr. (62) 

Here we introduced the notation (see Eq. (50)) 

Wo = — - — . (63) 

a 

Now, let t^ be the time moment when the contact force (49) reaches its maximum, i. e., 
P{i%) = 0. Then, by analogy with the modified incomplete loss angle Ss/2{^)^ we put 

4°/2H = f-^^M- (64) 

In view of (24), formula (64) can be rewritten in the form (30). 

The quantities £"3/2(01;, tug) and S^^2{^) (with explicit dependence on wq hidden) will be 
simply called the modified storage modulus and modified loss angle. 

It should be emphasized that the quantities E^^i^^^o) ^^^ ^3/2(^) depend on the layer 
thickness, although this fact is not refiected in the notation. Thus, in the viscoelastic 
case, the application of the indentation scaling factor /^s(a^) in the same way as on the 
right-hand side of formula (61) does not completely accounts for the thickness effect. 

Fig. 6a presents the comparison of the quantity E^^^i^^ ^0) with the modified incomplete 
storage modulus £3/2 (^) in the case of standard solid model. As it could be expected, 
the difference tends to zero as cj ^ 0, that is as the indentation protocol approaches the 
quasi-static limit. In the high-frequency range (as u) -^ oc), it can be also established 
rigorously that E^^^i^^^o) and Es/2{(^) both tend to £"(0) as well as Ei{(jj) does. On 
the contrary, the deference between S%{(jj) and 5m (^) does not vanish as cj ^ oc (see 
Fig. 6b), and the limit value of the relative error depends on the value of zuq. 
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Figure 6. Difference between E^^^{(jJ^vcjo)^ S^^^i^) ^^^ E^/2{(jj)^ ^3/2('^) fo^ different levels of 
indentation. 

3.2 Asymptotic analysis of E^^i^^^o) ^^ the low- and high-frequency limits 



To fix our ideas, let us assume that the relaxation modulus E{t) is determined by the 
Prony series 

t 



Hi / 

E{t) = E^ + J2Ejexp(- 



Pj 



(65) 



where Ej and pj are positive constants representing the relaxation strengths and relaxation 
times. Without loss of generality we may assume that Pi < P2 < • • • < Pm- 

Integrating by parts in (62) and changing the integration variable, we rewrite Eq. (62) as 



7r/2 

E^si^.wo) = E{0) + -- — r / E'(-]{co8zyf{woV^^^)dz. 



(66) 



Here again it is assumed that /3 = 3/2. 



In the low- frequency limit, the behavior oi E^^{uj^wq) as cj -^ will depend on the 
asymptotic behavior of the integral 



7r/2 

P^{uj,zuq) = j expf- 



ujpj 



(cos z)^ f{wo^/co8z) dz. 



(67) 



Using the formula (Dwight (1961), formula (567.9)) 

/ X exp(ax) ax = exp(ax) 2^ 
it can be easily shown that 



X 



k=0 



k\a^+^-^ ' 



7r/2 

/ 



exp 



UJpj J 



2ujp 



^n\{ujpjY^^-^ fiT\^ 



J^ k=0 



k\ 



2) 



= n!(cjp,)"+^ + 0('(a;p,)"+^ exp 



n 



2u>pj 



LJ ^0. 



(68) 



Thus, making use of formula (68) and the two-term Taylor expansion 

r(cosz)^/(^oVcosz) = \-z \^^^^ v-r\ +0(z ), 

one can now expand the integral /^(cj, tug) determined by (68) in power series with respect 
to up' . In such a way, we arrive at the following asymptotic representation: 



El{u^ ^0) = E^^u? £ E,p] (/3 + ^f ) + 0{u.% 



(69) 



Here it was taken into account that 



E{^) = E^^Y.^r 



(70) 



On the other hand, the following two-term asymptotic expansions holds true for the 
storage modulus (Tschoegl, 1997): 

Er{yj) = E^+uj^Y.EjP] + 0{uj^), uj^O. (71) 

Comparing asymptotic expansions (69) and (71), we see that the second terms on their 
right-hand sides coincide only if /5 = 1 and tuo = 0, when E^^{uj^W{)) coincides with the 
incomplete storage modulus Ei{(jj). 

In the high-frequency limit, the behavior of £"^(6^, tuo) as cj -^ oc depends on the smooth- 
ness properties of the function E'{t) at the point t = 0. In view of (65), Eq. (66) readily 
yields 

El[u, zuo) = EiO) + ^^f^ + Oi.-% (72) 

where E\0) = - E^^^ Ej/pj and 



7r/2 

r^(^o) = / {co8 z)^ f{zuo^/co8z) dz. 



On the other hand, the following two-term asymptotic expansions holds true for the 



storage modulus (Argatov, 2012): 



EAuj) = E(0) - ^^ + 0(cj-^), uj^ 00. 



UJ^ 



(73) 



Thus, based on the asymptotic expansions (72) and (73), it is estabhshed that E^{(jJ^zuo) 
tends to £"(0) as cj ^ oc as weU as Ei{(jj) does. 



3.3 Asymptotic analysis of 5^{(ju) in the low- and high-frequency limits 



Differentiating both sides of Eq. (49) with respect to time and taking into account the 
sinusoidal protocol (23), we get 



Ah'^Wn at ^ 



r d 
+ / — [{sin corf f{woV sin ujT)}E'{t - r) dr. 



(74) 



Substituting now the value (see Eq. (64)) 

i' 



4O _ ^ f^ ?0 



^M 



UJ\2 



into the equation P(i%) = in view of (74), we arrive at the foUowing equation: 



^/2-SO 



£;(0)a;sin^°(cos^°)'"V°(^°)= / ^{(cos(^? + 



xf{zuoyJcosCS$ + z))JE'(^-jdz. 
Here we introduced the notation 



(75) 



(76) 



Let C{6p^uj) and 7^.(52, cj) denote the left and right hand sides of Eq. (75), respectively. 



Following Argatov (2012), we construct solutions to Eq. (75), assuming that S^{u)) c^ CqU) 
as cj ^ and S^{(ju) c^ Coq/uj as cj ^ cx), where Co and Coo are constants. In both cases, 
^m(^) ^ 1 ^ViOh that 



C{5lu) = E{0)u~51:fI{Q) + 0[mf), ~5l ^ 0, 
where according to (76) we have 



(77) 



^0 ,1, 



J^(0)=/3/(t^o) + ^/Vo) 



In the low- frequency limit, making use of the asymptotic formula (68), we get 

m m 

J=l 3=1 

From (77) and (78), it follows that 

(jj *" 
5l{^)^T^Y.EjPj + 0{uj'), u^O. (79) 

Here the relation (70) was taken into account. 

We emphasize that the asymptotic representation (79) is in complete agreement with the 
leading term of the asymptotic expansion for the loss angle 5{(ju) as cj ^ 0. 

In the high-frequency limit, we will have 

7r/2 



/n 
— {{coszff{wo^c^sz)}E'{0)dz + 0{uj-^), w ^ oo. (80) 



Now, from (77) and (80), it follows that 



^>'-liw --■ <"' 



On the other hand, the following asymptotic representation holds true for the loss angle: 

Comparing relations (81) and (82), we see that they coincide only if /3 = 1 and f{^o) = 1. 
In the general case, in view of (76), we have 

Sjuj) - ^"(c.;) 2{/3 - l)f{TUo) + TUofizuo) 

Thus, according to (83), in the case of a viscoelastic half-space, when f{zuo) = 1, the 
relative error of the approximation Sp{(jj) for 5{(jj) approaches the value (/5 — l)//5 • 100%. 



4 Discussion 



The new material characteristics introduced above, that is the incomplete storage modulus 
Ei{(jj)^ the modified incomplete storage modulus Ei^{ijj)^ and the modified storage modulus 



Eq{(jj^wo)^ can be represented as follows: 



7t/{2uj) 



/d 
E{s)--{co8(jjs} 
as 





ds, 



7t/{2uj) 



M^) 



d 



E{s)^{{cosujsf]ds, 



7t/{2uj) 



E^^{u,wo) = - / E{s)-— I {cos ujs) 



J{woy/c08UJs)' 



ds. 



(84) 

(85) 
(86) 



In the same way, the storage modulus Ei{(jj) is recast as 



oo 

/d 
E{s)— -{cos (jjs} ds. 
ds 



(87) 



Thus, comparing formulas (84) - (86) with (87), we see that Ei{uj)^ Ef^{uj).^ and E^^{uj^ tug), 
represent a hierarchy of approximations for Eiiuj). Applying an asymptotic modeling 
approach for analyzing the interrelations between the new quantities, we have shown that 
the modified storage moduli asymptotically coincide with the storage modulus in the low- 
and high-frequency ranges. 

The values (3 = 1^ (3 = 3/2, and (3 = 2 correspond respectively to the cases of cylindri- 
cal, spherical, and conical indenters. The latter case also applies to pyramidal indenters 



(Giannakopoulos, 2006; Argatov, 2011). 



Fig. 7 illustrates the relationship between the parameters of viscoelastic materials mea- 
sured in a vibration indentation test and in a single indentation test with a fiat-ended 
cylindrical indenter. Due to the nonmonotonic behavior of the modified incomplete storage 
modulus Ei^{uj) with respect to the storage modulus Ei{uj) as it was shown in Fig. 4 (for 
the standard viscoelastic solid model) , the relationship between the parameters measured 
in the vibration and indentation tests with a spherical indenter will be more complicated. 

E2 t Measured in vibration test 




IVIeasured in indentation test 



Figure 7. Complex dynamic modulus schematic diagram. 



Now, let us consider the application of the developed theory to experimental data (Ronken 
et al., 2011). The experimental setup was described in detail elsewhere (Wirz et al., |2008 ). 
The modified moduli Ey2{^^ ^o) ^^d loss angles Sy2{^) of swine hyaline cartilage at two 
different locations calculated using the following formulas (see, Eqs. (61) and (30)): 



^3/2 (^^^o) = 



3(1 - jy') PjiJ 



7r{ir, 



Z 6^ 



M) 



(89) 



Here, /(tug) = /^s(<^m) is the indentation scaling factor corresponding to the maximum 
indentation depth and calculated according to Eq. (47). A Poisson's ratio of 0.5 was 
assumed. Note that the symbol t^ now denotes the time moment of maximum indentation 
instead of the symbol t^, because in the impact tests the indentation variation w{t) does 
not follow the sine law (23) precisely. 

Table 2 

Mean and standard deviation of the main parameters of two impact indentation tests (ten 
repetitions on one spot) for two swine cartilage samples with a spherical impactor of radius 
R = 0.5 mm and mass tti = 1.9 g (Ronken et al., 2011). 





Lateral condyle 


Medial condyle 


Sample thickness, h (nini) 


1.7 


1.9 


Initial indenter velocity, -uq (m/s) 


0.249 ± 0.002 


0.266 ± 0.002 


Time to maximum contact force, Im (ms) 


0.68 ± 0.01 


0.73 ± 0.01 


Maximum contact force, P{iM) (N) 


1.04 ±0.02 


0.96 ±0.01 


Indentation duration, im (ms) 


0.74 ±0.01 


0.81 ±0.01 


Maximum indentation, wq (mm) 


0.125 ±0.002 


0.145 ±0.003 


Contact force at maximum indentation, P{tm) (N) 


1.02 ±0.02 


0.93 ±0.01 


Effective angular frequency, uj (xlO^rad/s) 


2.13 ±0.02 


1.93 ±0.02 


Level of indentation, vuq 


0.147 ±0.001 


0.158 ±0.002 


Modified storage modulus, £'3/2(^,^^0) (MPa) 


15.3 ±0.5 


11.0 ±0.3 


Modified loss angle, 5^,^{(jo) (rad) 


0.114 ±0.007 


0.167 ±0.008 


Coefficient of restitution, e^, 


0.777 ±0.006 


0.722 ±0.012 



The data shown in the upper part of Table 2 was directly assessed in experiments, while 
the lower part of the table displays results evaluated according to the theory developed 
herein. The example illustrates the fact that the introduced characteristics £'3/2(6^,1170) 

and ^3/2(6^) depend on the effective angular frequency cj, which in turn depends on the 
initial indenter velocity v^ as well as on the mechanical properties of the sample itself. 
As it could be expected, at high frequencies, the modulus E^I2{uo^wq) increases with 

increasing cj, while the angle 5^i2{^) decreases (see Table 2). Note also that this example 

demonstrates a correlation in behavior of the modified loss angle 5^/2{^) a^d the coefficient 
of restitution e^. Of course, the impact indentation test requires a special consideration. 



but the observed characteristic behavior of £"3 /2(cl;, zuq) and S^^2{^) with the change in u) 
is quite typical. 

Now, let us discuss the significance of the developed mathematical approach from the 
viewpoint of formulating criteria for evaluation the quality of articular cartilage. In the 
cylindrical (fiat-ended) and spherical dynamic indentation tests, the following cartilage 
stiffness-related characteristics can be evaluated (see Eqs. (26) and (61), respectively): 

V ^^ (90) 



Here, Wq = w{tm) is the maximum indentation depth, P{tm) is the contact force corre- 
sponding to the time moment t = tm^ when the indenter reaches its maximum indentation 
depth. The criteria (90) and (91) in their static form (with no attention paid to the dy- 
namic nature of indentation process) have been used in a number of experimental studies 
on detection of degenerative changes in joint cartilage. 

First of all, it should be noted that the indentation scaling factors /^c(<^) and /^s(<^m) were 
evaluated under the assumption of isotropy and homogeneity of articular cartilage layer. 
It is anticipated that the inhomogeneity effect will be smaller in spherical indentation. 
In view of the layered structure of articular cartilage, the the anisotropy effect requires 
tacking into account the adjusted value of indentation scaling factor and the corresponding 
Poisson's ratio. Because Poisson's ratio u enters formulas (90) and (91) not only through 
the factor \ — v^ but also through the dependence of k^[ol) and /^s(<^m) on ^5 the question 
of the appropriate choice of v for the criteria (90) and (91) is more than academic, and it 
still remains open. 

Second, in dynamic indentation testing, the criteria (90) and (91) will depend on the 
loading protocol employed, because articular cartilage exhibits viscoelastic and poroelastic 
properties. Thus, in order to increase the sensitivity of the measurements with a hand- 
held indentation probe, the indentation protocol should be reproducible as well as the 
indentation time t^ should be kept the same. 

Third, in the framework of linear viscoelasticity, the criteria (90) and (91) are interpreted 
as the incomplete storage modulus E\{ijj) and the modified storage modulus £"3/2(6^,^170). 
Due to the linearity of the theory under consideration, the criterium (90) does not depend 
on the level of indentation. Hence, this fact could be used to check whether the linearity 
assumption is appropriate for small indentation depths, when w^jh < 0.1 or even less. On 
the other hand, the criterium (91) does depend on the level of indentation determined by 
the parameter (see Eq. (63)) 

Wq = — - — . 
h 

However, since the main manifestation of the thickness effect has been taken into account 

by means of the indentation scaling factor /^s(a^) evaluated at the maximum indentation 

depth, the dependence of E^I2{(jO^W{)) on w^) is rather weak as it is predicted by the 

standard viscoelastic solid model (see Fig. 6a). 



Further, in both indentation tests, one can also measure a dimensionless quantity that is 
directly related to time-dependent energy dissipation due to viscoelastic and poroelastic 
relaxation. Namely, the incomplete loss angle 6{(jj) and the modified loss angle ^3/2(^) 
were introduced based on the time shift between maxima of of the input (indentation dis- 
placement) and the output (contact force) through Eq. (30). It is important to emphasize 
that the quantity S{u))^ which is measured in the fiat-ended indentation test, does not 
depend on the thickness of the articular cartilage layer. At the same time, the thickness 
effect plays a crucial role in manifestation of the time-dependent response to indentation 
with a spherical indenter (see Fig. 6b). 

Finally, the developed viscoelastic models of cylindrical and spherical dynamic indentation 
tests allow one to compare the diagnostics criteria (that is diagnostics characteristics) 
experimentally measured by different indentation probes utilizing different indentation 
protocols (e.g., more closely approximating the actual movement of an operator's hand) 
as well as operating in different modes (vibration, dynamic indentation, impact testing). In 
view of the established fact that the incomplete storage modulus Ei{uj) and the modified 
storage modulus E^^^i'^i^o) asymptotically coincide with the storage modulus Ei{u)) in 
the low- and high-frequency ranges, the interrelationships between the different tests will 
be of particular interest in the middle frequency range. 



5 Conclusions 



We considered frictionless fiat-ended and spherical sinusoidally-driven indentation tests 
utilizing displacement-controlled loading protocol. In order to perform a rigorous analysis, 
we modeled the deformational behavior of articular cartilage tissue in the framework of 
viscoelasticity with a time-independent Poisson's ratio. In the linear case of fiat-ended 
indentation test, evaluating the dynamic indentation stiffness at the test turning point 
t = t^, we introduced the incomplete storage modulus Ei{(jj) for the effective frequency 
uj = 7T/{2tjn)- Considering the time difference t^ — Im between the time moments when 
the contact force reaches its maximum (dynamic stiffness vanishes at t = ^m) and the 
indenter displacement reaches its maximum (dynamic stiffness becomes infinite at t = t^), 
we introduced the so-called incomplete loss angle 5{uj). 

Analogous quantities were introduced in the nonlinear case of spherical sinusoidally-driven 
indentation test. First, when the sample thickness effect can be neglected, we introduced 
the modified incomplete storage modulus £"3/2 (^) and the modified incomplete loss angle 
^3/2 (^) (we use the same notation as in the linear case). Second, to account for the thick- 
ness effect, we introduced the indentation scaling factor /^s(<^) for the spherical indenter 
depending on Poisson's ratio and the relative contact radius a = a/h. Making use of the 
indentation scaling factor corresponding to the maximum indentation depth /^s(a^), we 
introduced the modified storage modulus E^^^i^^^o) which depends on the level of in- 
dentation characterized by the parameter wq = ^/w^/h. The modified loss angle S^^2{^) 
was introduced in the same way. 

We applied an asymptotic modeling approach for analyzing the interrelations between 
the new quantities E^^i^^ ^o)^ ^3/2 (^) ^^^ the classical characteristics Ei{uj)^ 6{(jj) in the 



low- and high-frequency ranges. It was shown that the modified storage modulus asymp- 
totically coincides with the storage modulus in the both limit cases, that is E^^i^^ ^o) — 
Ei{(jj) as cj ^ and uj ^ oo. However, coinciding with the loss angle 6{(jj) in the low- 
frequency range, the modified loss angle S^{(jj) markedly differs from S{(jj) in the high- 
frequency limit. We illustrated these facts for the standard viscoelastic solid model. 

Finally, the present study suggests that the use of dynamic indentation tests is largely 
twofold: the criteria (90) and (91) yield dimensional diagnostics characteristics, which 
could be related to some integral measure of material properties of the tested biological 
tissue; the second aspect of dynamic indentation diagnostics hinges on the importance 
of continuous monitoring of the tissue response to indentation. It is believed that both 
characteristics Ei{uj) (evaluated according to the criterium (90)) and S{(jj)^ which are as- 
sociated with fiat-ended indentation tests, can elicit perceptions of the mechanical quality 
of articular cartilage. In the dynamic non-destructive testing with a spherical indenter, in 
view of the fact that the thickness of a biological tissue sample is supposed to be unknown, 
the pronounced effect of the sample thickness on the modified loss angle S^^2{^) observed 
for different levels of indentation can be used as an indicator of the importance of the 
thickness effect for the modified storage modulus E^^i^^^o)- 



Acknowledgements 



One of the authors (LA.) gratefully acknowledges the support from the European Union 
Seventh Framework Programme under contract number PIIF-GA-2009-253055. 



References 

Appleyard, R.C., Swain, M.V., Khanna, S., Murrell, G.A.C., 2001. The accuracy and reli- 
ability of a novel handheld dynamic indentation probe for analysing articular cartilage. 
Phys. Med. Biol. 46, 541-550. 

Argatov, I.I., 2001. The pressure of a punch in the form of an elliptic paraboloid on an 
elastic layer of finite thickness. J. Appl. Math. Mech. 65, 495-508. 

Argatov, I.I., 2002. Characteristics of local compliance of an elastic body under a small 
punch indented into the plane part of its boundary. J. Appl. Mech. Techn. Phys. 43, 
147-153. 

Argatov, I.I., 2010. Frictionless and adhesive nanoindentation: Asymptotic modeling of 
size effects. Mech. Mater. 42, 807-815. 

Argatov, I., 2011. Depth-sensing indentation of a transversely isotropic elastic layer: 
Second-order asymptotic models for canonical indenters. Int. J. Solids Struct. 48, 3444- 
3452. 

Argatov, I., 2012. Sinusoidally-driven fiat-ended indentation of time-dependent materials: 
Asymptotic models for low and high rate loading. Mech. Mater. 48, 56-70. 

Argatov, I., Mishuris, G., 2011. An analytical solution for a linear viscoelastic layer loaded 
with a cylindrical punch: Evaluation of the rebound indentation test with application 
for assessing viability of articular cartilage. Mech. Res. Comm. 38, 565-568. 



Armstrong, C.G., Lai, W.M., Mow, V.C, 1984. An analysis of the unconfined compression 
of articular cartilage. J. Biomech. Eng. 106, 165-173. 

Bae, W.C, Temple, M.M., Amiel, D., Coutts, R.D., Niederauer, G.G., Sah, R.L., 2003. 
Indentation testing of human cartilage: Sensitivity to articular surface degeneration. 
Arthritis Rheum. 48, 3382-3394. 

Borodich, F.M., Keer, L.M., 2004. Contact problems and depth-sensing nanoindentation 
for frictionless and frictional boundary conditions. Int. J. Solids Struct. 41, 2479-2499. 

Broom, N.D., Flachsmann, R., 2003. Physical indicators of cartilage health: The relevance 
of compliance, thickness, swelling and fibrillar texture. J. Anat. 202, 481-94. 

Cao, Y., Ma, D., Raabe, D., 2009. The use of fiat punch indentation to determine the 
viscoelastic properties in the time and frequency domains of a soft layer bonded to a 
rigid substrate. Acta Biomaterialia 5, 240-248. 

Cheng, Y.-T., Yang, F., 2009. Obtaining shear relaxation modulus and creep compli- 
ance of linear viscoelastic materials from instrumented indentation using axisymmetric 
indenters of power-law profiles. J. Mater. Res. 24, 3013-3017. 

Christensen, R.M., 1971, Theory of Viscoelasticity. Academic Press, New York. 

de Freitas, P., Wirz, D., Stolz, M., Gopfert, B., Friederich, N.-F., Daniels, A.U., 2010. 
Pulsatile dynamic stiffness of cartilage-like materials and use of agarose gels to validate 
mechanical methods and models. J. Biomed. Mater. Res. Part B: Appl. Biomater. 78B, 
347-357. 

Dwight, H.B., 1961. Tables of Integrals and Other Mathematical Data. The Macmillan 
Company, New York. 

Galin, L.A., 1946. Spatial contact problems of the theory of elasticity for punches of 
circular shape in planar projection. J. Appl. Math. Mech. (PMM) 10, 425-448 (in 
Russian) . 

Giannakopoulos, A.E., 2006. Elastic and viscoelastic indentation of fiat surfaces by pyra- 
mid indentors. J. Mech. Phys. Solids 54, 1305-1332. 

Hayes, W.C, Keer, L.M., Herrmann, G., Mockros, L.F., 1972. A mathematical analysis 
for indentation tests of articular cartilage. J. Biomech. 5, 541-551. 

Hayes, W.C, Mockros, L.F., 1971. Viscoelastic properties of human articular cartilage. J. 
Appl. Physiol. 31, 562-568. 

Hu, K., Radhakrishnan, P., Patel, R.V., Mao, J.J., 2001. Regional structural and vis- 
coelastic properties of fibrocartilage upon dynamic nanoindentation of the articular 
condyle. J. Struct. Biol. 136, 46-52. 

Korhonen, R.K., Saarakkala, S., Toyras, J., Laasanen, M.S., Kiviranta, I., Jurvelin, J.S., 
2003. Experimental and numerical validation for the novel configuration of an arthro- 
scopic indentation instrument. Phys. Med. Biol. 48, 1565-1576. 

Kren, A. P., Naumov, A.O., 2010. Determination of the relaxation function for viscoelastic 
materials at low velocity impact. Int. J. Impact Eng. 37, 170-176. 

Kusano, T., Jakob, R.P., Gautier, E., Magnussen, R.A., Hoogewoud, H., Jacobi, M., 2011. 
Treatment of isolated chondral and osteochondral defects in the knee by autologous 
matrix-induced chondrogenesis (AMIC). Knee Surg. Sports Traumatol. Arthrosc. DOI 
10.1007/s00167-011-1840-2. 

Lau, A., Oyen, M.L., Kent, R.W., Murakami, D., Torigaki, T., 2008. Indentation stiffness 
of aging human costal cartilage. Acta Biomater. 4, 97-103. 

Lyyra-Laitinen, T., Niinimaki, M., Toyras, J., Lindgren, R., Kiviranta, I., Jurvelin, J.S., 
1999. Optimization of the arthroscopic indentation instrument for the measurement of 
thin cartilage stiffness. Phys. Med. Biol. 44, 2511-2524. 



Oyen, M.L., 2005. Spherical indentation creep following ramp loading. J. Mater. Res. 20, 
2094-2100. 

Parsons, J.R. and Black, J., 1977. The viscoelastic shear behavior of normal rabbit artic- 
ular cartilage. J. Biomech. 10, 21-29. 

Pipkin, A.C., 1972. Lectures on Viscoelasticity Theory. Springer, Berlin. 

Ronken, S., Arnold, M.P., Ardura Garcia, H., Jeger, A., Daniels, A.U., Wirz, D., 2011. 
A comparison of healthy human and swine articular cartilage dynamic indentation 
mechanics. Biomech. Model. Mechanobiol. DOI: 10. 1007/sl0237-01 1-0338-7. 

Schinagl, R.M., Gurskis, D., Chen, A.C., Sah, R.L., 1997. Depth-dependent confined com- 
pression modulus of full-thickness bovine articular cartilage. J. Orthop. Res. 15, 499- 
506. 

Suh, J.-K., Li, Z., Woo, S.L.-Y., 1995. Dynamic behavior of a biphasic cartilage model 
under cyclic compressive loading. J. Biomech. 28, 357-364. 

Ting, T.C.T., 1968. Contact problems in the linear theory of viscoelasticity. J. Appl. 
Mech. 35, 248-254. 

Tschoegl, N.W. 1997. Time dependence in material properties: An overview. Mech. Time- 
Depend. Mat. 1, 3-31. 

Toyras, J., Lyyra-Laitinen, T., Niinimaki, M., Lindgren, R., Nieminen, M.T., Kiviranta, 
I., Jurvelin, J.S., 2001. Estimation of the Young's modulus of articular cartilage using 
an arthroscopic indentation instrument and ultrasonic measurement of tissue thickness. 
J. Biomech. 34, 251-256. 

Vorovich, LL, Aleksandrov, V.M., Babeshko, V.A. 1974. Non-classical Mixed Problems of 
the Theory of Elasticity. Nauka, Moscow [in Russian] . 

Wirz, D., Kohler, C, Keller, K., Gopfert, B., Hudetz, D., Daniels, A.U., 2008. Dynamic 
stiffness of articular cartilage by single impact micro-indentation (SIMI). J. Biomech. 
41, Suppl. 1, P. S172. 

Zhang, C.Y., Zhang, Y.W., 2004. Extracting the mechanical properties of a viscoelastic 
polymeric film on a hard elastic substrate. J. Mater. Res. 14, 3053-3061. 



